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Abstract 

The effect of pitch-angle anisotropy of fast electrons on generation of nonthermal radio emission is stud- 
ied. Incoherent gyrosynchrotron radiation is shown to depend strongly on the anisotropy. In particular, 
the spectral index of gyrosynchrotron radiation increases up to a factor of 3-4 compared with the isotropic 
case. The degree of polarization (X-mode) increases for loss-cone distributions, while decreases for beam- 
like distributions. Ample evidence of the pitch- angle anisotropy effect on (1) spatial distribution of the 
radio brightness, (2) spatially resolved light curves of the intensity and polarization, (3) spectral hardness 
of microwave bursts, is found by exploring observations performed with Nobeyama Radioheliograph and 
Owens Valley Solar Array. Above some threshold in the angular gradient, the electron cyclotron maser 
instability (coherent gyrosynchrotron emission) can develop, provided that standard gyrosynchrotron emis- 
sion is accompanied by a lower-frequency intense coherent emission. The coherent emission is studied in 
detail for a realistic distribution of fast electrons over the energy and pitch-angle. In agreement with 
observations of narrowband radio spikes, it is found that (1) hard (power-law) distributions over energy 
are preferable to produce the coherent emission, (2) the threshold of the instability corresponds to quite 
an anisotropic electron distribution, thus, the pitch- angle anisotropy derived from the properties of the 
continuum gyrosynchrotron radiation will not necessarily give rise to coherent emission (either enhanced 
isotropization implied by quasilnear saturation of the instability). 

Key words: acceleration of particles — instabilities — radiation mechanisms: non-thermal — Sumflares — 
Sumradio radiation 



1. Introduction 

Observations of solar radio bursts with Nobeyama 
Radioheliograph (NoRH) revealed unambiguously the 
presence of anisotropic pitch-angle distributions of fast 
electrons accelerated in flares and accumulated in the ra- 
dio sources. In particular, the existence of the microwave 
loop-top brightness peak of optically thin gyrosynchrotron 
(GS) emission (Melnikov et al. 2001; Kundu et al. 2001) 
has been interpreted as tracer of a strong concentration of 
mildly relativistic electrons trapped and accumulated in 
the top of a flaring loop (Melnikov et al. 2002a), which is 
an evidence of the transverse anisotropy of radiating elec- 
trons. Then, Lee & Gary (2000) studied the evolution 
of the radio spectrum of 1993 Jun 3 flare observed with 
Owens Valley Solar Array (OVSA) and found important 
indication of the anisotropic injection of fast electrons into 
the radio source and their consequent isotropisation due 
to angular scattering. The importance of the pitch-angle 
anisotropy for the modern radio data interpretation has 
been widely discussed during this Nobeyama Symposium 
(Altyntsev 2005; Bastian 2005; Gary 2005; Melnikov 
2005). This paper discusses the effect of the pitch-angle 
anisotropy on generation of GS radiation (both incoherent 
and coherent). 

A detailed study (stimulated by the mentioned findings) 
of incoherent GS emission produced by anisotropic pitch- 
angle distributions of the loss-cone and beam-like types 
is given in recent papers (Fleishman & Melnikov 2003a; 



Fleishman & Melnikov 2003b). Accordingly, we outline 
only briefly why the anisotropy has actually a strong effect 
on GS emission. 

The effect of anisotropy on the optically thin emission 
is primarily related to the energy-dependent directivity of 
the GS radiation. Indeed, most of the synchrotron radia- 
tion by a single ultrarelativistic electron is emitted within 
a narrow cone 

d ~ 7- 1 = mc 2 / E (1) 

along the particle velocity, where 7 e is the Lorentz-factor 
of the electron, m, E are the mass and energy of the elec- 
tron, c is the speed of light, so {) <C 1 if -f e ^> 1. On the 
other hand, for semirelativistic and mildly relativistic elec- 
trons responsible for the GS emission, the directivity is not 
quite as strong, so the anisotropy (if present) can signifi- 
cantly change the radiation spectrum when the transition 
from moderately relativistic to highly relativistic regime 
occurs. Typical energies of fast electrons producing so- 
lar microwave continuum bursts (which are commonly be- 
lieved to be produced by GS mechanism, e.g., Bastian et 
al. 1998) range from hundreds keV to a few MeV. Thus, 
the pitch-angle anisotropy of fast electrons in solar flares 
is particularly important for generation of microwave con- 
tinuum radio emission, as well as for driving instabilities, 
which can result in coherent radio emissions. 
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Fig. 1. Top: Angular distribution functions: gaussian 
(Fleishman & Yastrebov 1994a), dotted line; sin-N loss-cone 
(Aschwandcn 1990) for N = 6, dash-dotted line; exponen- 
tial loss-cone (5) for A/i = 0.15, dashed line; and gaussian 
loss-cone (6) for A/i = 0.15, solid line. The distributions are 
plotted for no = cos0 c = 60°, all of them are normalized to 
unity. Middle: The derivatives of the angular distributions 
over fA. Bottom: the ratios of the derivatives to the distribu- 
tion functions. 

2. Angular Distribution Function 

As a first approximation we consider a (simplified) fac 
torized distribution function of fast electrons 
N, 
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The distribution function (2) is normalized by d 3 p 

J f(p)p 2 dpdLLd<p = N e , (3) 

where N e is the number density of fast electrons, /i = cosO, 
6 is the pitch-angle, ip is the azimuth angle, the specific 
distributions fi(p) and /2O?) are normalized as follows: 
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We assumed the distribution function over the momentum 
modulus to be a kind of power-law, while we consider a few 
different forms for the pitch-angle distribution to study 
the sensitivity of the GS properties to the shape of the 
angular distribution. 

In addition to the gaussian and sin-N functions (ex- 
plored by Fleishman & Melnikov 2003b), we use also the 
exponential loss-cone function 
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and the gaussian loss-cone function 
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Let us consider these angular distributions (presented 
in figure 1) in more detail, having in mind that along with 
modification of incoherent GS emission the anisotropy can 
give rise to some instabilities (which depends on relative 
value of the angular gradient), in particular, to electron 
cyclotron mascr (ECM) emission at the low gyroharmon- 
ics. The gaussian angular distribution provides rapid de- 
crease of the number of electrons toward the direction 
of the magnetic field, and affect the GS substantially 
(Fleishman & Melnikov 2003b). 

However, this distribution is not favorable to drive the 
instabilities (Fleishman & Yastrebov 1994a, 1994b), be- 
cause the instability related to the gaussian angular distri- 
bution is rather weak being only slightly above the thresh- 
old. The efficiency of the instability is set up by relative 
importance of negative derivative of the distribution func- 
tion over the momentum modulus and positive contribu- 
tion related to angular gradient, which can be quantita- 
tively described by the ratio of the angular gradient of the 
pitch-angle distribution to the angular distribution func- 
tion itself. Although the angular derivative in figure 1 is 
not small, the peak value occurs at those angles where the 
angular distribution function is large, the ratio of these 
values is, therefore, small everywhere, figure 1 bottom. 
Thus, the gaussian pitch- angle distribution can be stable 
against instabilities for rather broad range of parameters. 

Sin-N angular distribution (Aschwandcn 1990) has very 
different properties. While this function looks rather 
smooth (it is continuous function together with its first 
derivative) and, thus, attractive for the use in theoretical 
calculations, it cannot be considered as a good approx- 
imation for the angular distribution since this function 
has a very deep intrinsic problem. Indeed, if we compare 
the behavior of the function and its first derivative, we 
discover that the derivative decreases (with the decrease 
of the pitch-angle) much more slowly than the function 
itself. Particularly, the arrows in figure 1 indicate the re- 
gion where the derivative is of order of unity, while the 
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function is much less than unity. Note, that the ratio di- 
verges as the pitch-angle decreases. Thus, sin-N function 
describes the case of highly unstable electron distribution, 
which gives rise to ECM instability operating significantly 
above the respective threshold for any parameters of this 
function. As a result, the instability exists for almost 
any shape of the momentum distribution of fast electrons, 
and its efficiency has extremely weak dependence on the 
electron spectral index. Such highly unstable distribution 
function can hardly be built up under reasonable assump- 
tions about fast electron acceleration and transport. 

This artificial property has already produced a few un- 
expected and poorly understood results. In particular, 
the detailed numerical calculations of the ECM quasilin- 
ear saturation (Aschwanden 1990) did not ensure the en- 
tire maser saturation even after 50 growth times: the peak 
value of the growth rate was reduced by less than one order 
of magnitude for this time. The reason for this surpris- 
ing behavior is rather slow angular diffusion in the range 
of small pitch-angles, where the distribution function is 
small but the angular gradient is sufficiently large to pro- 
vide the instability and then support it for quite a long 
time. The quasilinear relaxation occurs differently if the 
instability operates closer to its threshold (Fleishman & 
Arzner 2000). 

Exponential loss-cone function (5) is free from such kind 
of problem. From the theoretical point of view it has only 
one disadvantage - discontinuity of its first derivative at 
6 = 9 C . Alternatively, a similar loss-cone distribution but 
with faster, gaussian, decrease of the number of electrons 
in the loss-cone (6) is also attractive from the theoretical 
point of view, since it is a function continuous together 
with its first derivative. 

3. Incoherent Gyrosynchrotron Emission 

The study of the GS emission produced by anisotropic 
pitch-angle distributions requires exact calculation of 
the corresponding emission and absorption coefficients 
(Fleishman & Mclnikov 2003b). This is especially im- 
portant for the directions along which a relatively small 
amount of particles is moving, so adding up the contri- 
butions from numerous electrons moving at unfavorable 
directions should be done correctly. Figure 2 presents an 
example of intensity, degree of polarization and spectral 
index of GS radiation produced by a gaussian angular dis- 
tribution vs frequency as the degree of the anisotropy (an- 
gular gradient) changes (Fleishman & Melnikov 2003b). 

The effect of the anisotropy on the properties of GS 
emission depends evidently on the viewing angle between 
the directions of wave propagation and the magnetic field. 
For the quasitransverse (QT) case, r\ = 0.2 (z? = 78°), 
right column, the effect of the pitch-angle anisotropy on 
the spectral index (bottom panel) is rather weak, while 
the intensity and degree of polarization increase with 
anisotropy increases. However, the anisotropy has exceed- 
ingly large effect on GS radiation produced at a quasi- 
parralel (QP) direction, r} = 0.8 (i9 = 37°), left column. 
Indeed, the radiation intensity changes by orders of mag- 
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Fig. 2. GS radiation intensity, degree of polarization, and 
spectral index vs. frequency for various /^o in the Gaussian 
pitch-angle distribution fa oc exp |^ — (/^ — Aii) 2 //^^] °f * nc 
loss-cone type, fi± = 

nitudc compared with the isotropic case. Respectively, the 
difference in the spectral index in the optically thin fre- 
quency range is up to a factor of four (left panel, bottom). 
Remarkable change of polarization occurs as well: the de- 
gree of polarization increases in the optically thin region 
and can approach 100% for sufficiently high anisotropy. 
Furthermore, the sense of polarization can correspond 
to X-modc in the optically thick region contrary to the 
isotropic case. 

Complementary, figure 3 shows the intensity, degree of 
polarization, and spectral index of GS emission observed 
at four different viewing angles from the gaussian (6) and 
exponential (5) loss-cone distributions (all other condi- 
tions are the same). GS emission is evidently very sen- 
sitive to the type of the angular distribution of fast elec- 
trons. The difference between these two distributions are 
seen most prominently in the spectral index plots: al- 
though the spectral index can be as large as 6 for the 
gaussian loss-cone distribution, it does not exceed 3 for 
the exponential loss-cone distribution. 

Thus, this high sensitivity of the results to the type 
of the angular distribution can, in principle, be used to 
obtain important constraints on the pitch-angle distribu- 
tions present during solar flares. We emphasize, that along 
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Fig. 3. GS radiation produced by fast electrons with a 
power-law distribution (with the index 7 = 7) over momen- 
tum modulus and gaussian (6) (left) or exponential (5) (right) 
loss-cone functions. The total intensity, polarization and local 
spectral index are shown in the plot. 

with strong effect on radiation intensity and spectral in- 
dex, the anisotropy influences the degree of polarization, 
which also might be applied for quantitative diagnostics 
of the pitch-angle distribution. 

Finally, we have to mention strong effect of pitch-angle 
anisotropy on the low- frequency "harmonic" structure of 
GS radiation (Fleishman & Melnikov 2003a). However, 
detection of this structure in the typical non-homogeneous 
sources requires high-resolution imaging spectroscopy, 
which is unavailable yet. 
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Fig. 4. Example of light curves of the microwave emission 
from a few positions along the loop for the case of isotropic 
injection function and angular-dependent life-time of fast elec- 
trons (top). The thin line represents the injection profile. 
Corresponding distributions (at a few time frames) of radio 
brightness (middle) and degree of polarization (bottom) along 
the magnetic trap. 



4. Anisotropy Build up in the Magnetic Loops 

The importance of transport effects (Kennel & 
Petscheck 1966; Melrose & Brown 1976) for genera- 
tion of nonthermal emissions during solar flares has been 
recognized long ago, in particular, the magnetic mirror- 
ing was argued to be a key process to account for the 
loop-top brightness peak of the optically thin GS emis- 
sion (Melnikov et al. 2002a). This finding implies that the 
anisotropy can actually be built up due to transport effects 
even if the injection of fast electrons occurs isotropically 
(although this does not deny the possibility of anisotropic 
injection of fast electrons). 

Typical treatment of the electron trapping at the mag- 
netic loops frequently explores the idea of empty loss-cone 
due to particle precipitation and losses at the dense foot- 
points in the weak diffusion limit. Although this idea 



brings an apparent advantage of the simplicity, it has also 
an important disadvantage of the discontinuity of the dis- 
tribution function over the pitch-angle resulting in steep 
angular gradient, which affects strongly the calculation of 
the absorption coefficients. 

To avoid the problems with the discontinuity and ac- 
count for magnetic mirroring, energy losses and weak an- 
gular diffusion, we put forward a new model solution of 
the transport equation: 

f(E,n,s,t) = f expf *~* )g(E,n,s,t')dt'{7) 

J-oo \ t{E,h,s)J 

where g(E, fj,,s,t') and f(E, fJ>,s,t) are the injection and 
distribution functions of the fast electrons at the radio 
source vs energy E, the cosine of pitch angle ^, position 
along the loop s, and time t . Here the phcnomcnological 
life-time parameter r can depend on the energy, position, 
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and pitch-angle. 

For the weak diffusion regime we adopt that the life 
time r display more or less fast decrease inside the loss- 
cone (when [i approaches the unity). This ensures smooth 
behavior of the distribution function, although the an- 
gular gradients can still be large for a favorable combi- 
nation of parameters, and the net anisotropy can grow 
with time. Examples of the model radio brightness dis- 
tributions at the subsequent times as a radio burst devel- 
ops, and the spatially resolved light curves (at an opti- 
cally thin frequency) in the context of anisotropic trans- 
port are shown in figure 4. To reveal the net effect of 
the angular-dependent transport, figure 4 assumes the life 
time of the electrons to depend on the pitch-angle of the 
electron only. It might be noted that the radio bright- 
ness peak moves towards the loop-top as the burst de- 
velops in agreement with observations (Melnikov et al. 
2002a). Complementary the emission from different lo- 
cations along the loop displays systematic delay increas- 
ing from the foot-points to the loop-top, which has also 
been widely observed with NoRH data (Melnikov et al. 
2002b; Melnikov 2005; Bastian 2005). 

The dependence of the degree of polarization on time 
and its spatial distribution arc highly affected by the 
pitch-angle anisotropy. Indeed, with isotropic injection, 
magnetic mirroring, and weak angular diffusion adopted 
in figure 4, the anisotropy increases with time progres- 
sively, which gives rise to the increase of the degree of 
polarization with time in the legs and top of the loop. At 
the foot-points (r/ro > 1.4), however, the angular distri- 
bution repeats the injected distribution since the life time 
of the electrons at the foot-points is too small to build up 
noticeable pitch-angle anisotropy. 

5. Observational Evidence of Anisotropic Pitch- 
Angle Distribution in Flares 

Besides already discussed spatial distributions of the 
radio brightness and spatially localized light curves, be- 
havior of the degree of polarization typical for anisotropic 
pitch-angle distributions is frequently observed with 
NoRH. An example of the localized light curve (from one 
pixel of the magnetic loop leg) and corresponding time 
profile of the degree of polarization for the flare observed 
on 28-May-1993 is shown in figure 5. Remarkably, the 
absolute value of the degree of polarization keeps grow- 
ing during the main peak of the burst, which is indicative 
for the increase of the anisotropy probably due to parti- 
cle losses through the loss-cone, since other parameters 
affecting the polarization (magnetic field and the viewing 
angle) cannot change much at the localized region of the 
loop. The increase of the degree of polarization is followed 
by its decrease at a later stage of the burst. Interestingly, 
this turning point in the degree of polarization plot occurs 
at the same time as the break in the light curve where 
rapid decay of the radiation intensity slows down notice- 
ably. This indicates that the losses due to precipitation 
into the loss-cone (which provides the anisotropisation of 
the electron distribution) are giving way to another decay 
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Fig. 5. A localized light curve (top) and and degree of po- 
larization (bottom) for the 28-May-1993 event as observed at 
17 GHz with NoRH (courtesy by T.S. Bastian). 

process, probably, to Coulomb losses of already accumu- 
lated fast electrons. The Coulomb losses result it harder 
electron spectra, which generate GS radiation with weaker 
polarization. 

A prominent spectral signature of the pitch-angle 
anisotropy has been found from the detailed analysis of 
two-component radio bursts composed of a high-frequency 
GS component and lower-frequency coherent component 
(Fleishman et al. 2003) observed with OVSA. Indeed, the 
microwave sources (accompanying the coherent dm spiky 
emission) observed at QP viewing angles to the magnetic 
field, figure 6, displays remarkably softer spectra at the 
optically thin region, than all other events on average, 
which can only be understood within the effect of loss- 
cone angular anisotropy on the GS spectrum. Moreover, 
all of the microwave bursts were highly polarized in the 
optically thin range (Fleishman et al. 2003), which is most 
probably related to the pitch-angle anisotropy of the loss- 
cone type (see figures 2, 3). 
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Fig. 6. The histograms of the high-frequency microwave 
spectral indices (Fleishman et al. 2003) for all bursts recorded 
by OVSA during 2001 (thin line), superimposed on the cor- 
responding histogram for spike-producing microwave bursts 
observed at QT viewing angle (left) and QP viewing angle 
(right), thick lines. Strong excess of the steep spectra is ob- 
served in the QP case. 

The difference in the spectral indices for these two 
classes of spike-producing microwave bursts (QT and QP) 
is as high as 3-4, which implies some important constraints 
on the shape of the angular distribution function. For ex- 
ample, the exponential decrease of the number of electrons 
in the loss-cone (5) turns to be too slow to provide the ob- 
served difference of 3-4 in spectral indices of optically thin 
GS emission at QP and QT directions (figure 3, right col- 
umn). However, the gaussian loss-cone distribution (6) 
can evidently provide the required difference in the opti- 
cally thin spectral indices (figure 3, left column). One can 
also note, that the gaussian loss-cone distribution provides 
stronger degree of polarization in the optically thin range 
than the exponential loss-cone distribution. 

6. Electron Cyclotron Maser Emission 

It is clear and has been well-known for a long time, 
that the pitch-angle anisotropy can drive instabilities, in 
particular, ECM instability. In essence, ECM emission is 
a coherent GS emission arising around lowest harmonics of 
the gyrofrcquency when the absorption coefficient changes 
the sign and turns to be the amplification coefficient. 

The theory of ECM emission in the solar corona has 
been extensively developed for the recent years (see 
Fleishman (2004a) for more detailed discussion and refer- 
ences). Big efforts have been made to consider the linear 
stage of the instability to specify the fastest growing wave- 
mode as a function of the plasma parameter Y — Lu pe /ujBe, 
where uj pe is the electron plasma frequency and tOBe is the 
electron gyrofrequency. 

Then, a particular attention has been paid to the role 
of absorption by the background coronal plasma in either 
emission source or along the ray path (in the harmonic gy- 
rolayers) (Sharma et al. 1982; Sharma & Vlahos 1984). 
The absorption in gyrolayers was found to represent a se- 



vere problem for the ECM emission to escape the corona 
(Sharma & Vlahos 1984), since the optical depth of the 
second gyrolayer calculated with the maxwcllian back- 
ground distribution is typically much larger than unity. 
However, estimates of nonlinear absorption (when the 
background distribution is modified by the incident pow- 
erful radiation) (Huang 1987 ) and particle-in-cell simula- 
tions (McKeen et al. 1990) show that even fundamental 
ECM emission has a chance to escape. Nevertheless, it 
remains unclear if the nonlinear absorption actually hap- 
pens in solar conditions or not (and, accordingly, which 
harmonics can really escape the source to the observer). 

Nonlinear stage of the ECM instability has also been 
studied. One-dimensional particle-in-cell simulations 
(Pritchett 1986) were found to agree well with a sim- 
plified analytical approach (Wu et al. 1981). The one- 
dimensional quasilincar saturation of the instability gives 
rise to a plateau formation in the distribution function of 
fast electrons over the transverse (to the magnetic field) 
velocity, providing, therefore, an absolute stabilization of 
the electron distribution. This means, that only the wave- 
mode with the largest growth rate can be amplified sub- 
stantially, while weaker growing wave-modes remain prac- 
tically unamplificd. 

The quasilincar relaxation is well known to occur rather 
differently for the two-dimensional case (Akhiezer et al. 
1974) relevant to the ECM emission. Indeed, the relax- 
ation due to the dominant mode affects the electron dis- 
tribution in such a way that the new (modified) distribu- 
tion is stable against this particular wave-mode. However, 
this new state is not necessarily stable against generation 
of other wave- modes, i.e., it is not absolutely stable dis- 
tribution. Thus, other wave-modes can well continue to 
grow until they start to affect the electron distribution in 
their turn. 

Two-dimensional fully relativistic numerical calcula- 
tions of the ECM quasilincar relaxation (Aschwanden 
1990) revealed many important properties of the relax- 
ation process against the fastest growing mode. However, 
the initial sin-N distribution function adopted for the 
calculations has artificially large angular gradient in the 
range of pitch-angles practically free from fast electrons 
as has been already discussed in figure 1. Quasilinear 
relaxation has further been studied numerically for a par- 
ticular set of the involved parameters when the lower- 
hybrid wave-mode has the largest growth rate (Fleishman 
& Arzner 2000). The electron distribution was shown 
to remain clearly anisotropic after the relaxation against 
the lower-hybrid waves. The new anisotropic distribu- 
tion is, however, stable against the lower-hybrid waves, al- 
though unstable against the electromagnetic wave-modes 
in agreement with the general results of the plasma physics 
(Akhiezer et al. 1974). In addition, the time profile of 
the generated wave energy density was found to resemble 
the observed spike time profiles (Giidel & Bcnz 1990). 
Nevertheless, the ECM saturation process has not been 
studied in detail yet, thus, additional studies of the quasi- 
lincar relaxation against electromagnetic wave-modes are 
necessary. 
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The role of coronal inhomogeneities both in ECM source 
and at the ray path was analyzed. Effect of MHD- 
waves on the longitudinal transparency window (Robinson 
f991a) and on the efficiency of the mode conversion 
(Robinson 199fb) was studied. Then, ECM line forma- 
tion in a source with random inhomogeneities was found to 
produce either spectral broadening or splitting of the ini- 
tial peak depending on the strength of the inhomogeneity 
(Fleishman 2004b) . Mutual action of regular coronal non- 
uniformity and random local inhomogeneities was consid- 
ered for a beam-like electron distributions (Vlasov et al. 
2002). 

In the previous section we concluded that the gaussian 
loss-cone function (6) is a good candidate for the electron 
pitch-angle distribution in the main source of continuum 
GS emission in the spike-producing events. Let us con- 
sider here the main properties of ECM emission produced 
by fast electrons with this angular distribution and power- 
law distribution over momentum modulus to be compared 
with the observed properties of the radio spikes appeared 
at the main phase of the burst. The results presented 
below are obtained numerically from fully relativistic cal- 
culations of the exact expressions of the GS absorption 
coefficient. 

6.1. Basic parameters. 

The growth rate value depends on many parameters, 
therefore, we first selected a set of basic parameters and 
then varied each of them successively. This approach al- 
lows to discuss properly what effect is produced by each 
particular parameter. 

For our calculations we selected two basic values of the 
plasma frequency to gyrofrcquency ratio, Y = Lo pe /ujBe = 
1, 1.3. When Y = 1, the most rapid growth of the ex- 
traordinary waves occurs at the second harmonics, s = 
u>/u>Be ~ 2 and of the ordinary waves at the fundamental, 
s«l, while for Y = 1.3 the growth of 01 waves becomes 
ineffective, and it is giving way to 02 waves. 

The angular distribution (6) depends on two param- 
eters: loss-cone angle 8 C and the width of the distribu- 
tion Afi. We selected 8 C = 60° and = 0.15. The later 
value allows the instability to operate in the moderately- 
above-threshold regime. For the basic distribution of fast 
electrons over momentum we selected the following pa- 
rameters typical for flares: po/mc = 0.2 (E^m ~ 10 keV), 
Pb r /mc = 3 (Ekin ~ 1.1 MeV) and 7 = 6. It is assumed 
that at p > pi, r the power-law distribution is giving way 
to exponential cut-off. 

6.2. Effect of plasma density. 

The effect of plasma density on the growth rates at 
the lowest harmonics of the gyrofrequency was repeat- 
edly studied for various electron distributions (Sharma et 
al. 1982; Sharma & Vlahos 1984; Li 1986; Winglee & 
Dulk 1986; Aschwandcn 1990; Fleishman & Yastrebov 
1994a; Fleishman & Yastrebov 1994b; Ledcncv 1998; 
Stupp 2000; Fleishman 2004a). It has been established 
that for small Y < 0.2 — 0.4 the fundamental extraordinary 
waves (XI) have the largest growth rates. Then, for larger 
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Fig. 7. Dimcnsionlcss spatial growth rates (normalized by 
^ 2c " -jf- ) for three first harmonics of the gyrofrequency for 
various parameters: 7 = 6 (left column), 7 = 8 (right column); 
e c = 45° (top), 6c = 60° (middle), 6 C = 75° (bottom). Set of 
basic parameters is used: A^t = 0.15, po/mc = 0.2, p br /mc = 3. 



Y , the amplification of XI waves becomes inefficient and 
other waves become dominant (01, X2, ...). Our calcu- 
lations agree well with this standard picture (since it is 
largely related to kinematic restrictions) . 

The dependence of the growth rates on Y (for the 
adopted power-law momentum distribution and gaussian 
loss-cone angular distribution) in the limits 0.1 < Y < 2.5 
is given in figure 7. The growth rate of each harmon- 
ics decreases rapidly when the cut-off frequency of cor- 
responding wave-mode exceeds the respective multiple of 
gyrofrequency (i.e., when Y > 1 for 01 waves, Y > V2 for 
X2 waves, Y > 2 for 02 waves etc). The absolute peak of 
the growth rate at a particular harmonics of the extraor- 
dinary waves is noticeably larger than of ordinary waves 
at the same harmonics. The absolute peak value decreases 
as the harmonic number increases. 

The range of Y variation where a wave-mode is the 
fastest growing mode depends on parameters of the fast 
electron distribution. In particular, for small loss-cone an- 
gle value (top), the growth rates of 01 waves are larger 
than of X2 waves for Y < 1. However, the difference be- 
tween them decreases as 9 C increases (middle) and X2 
waves grow faster than 01 waves for large 8 C (bottom) 
in agreement with the finding by Winglee & Dulk (1986). 
Respectively, 01 mode grows faster than X2 mode if Y < 1 
for 9 C = 45°, if Y < 0.85 for C = 60°, and never for 6 C = 75°. 
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Fig. 8. Spatial growth rates for X2-, Ol-, and 02- modes 
vs width of the loss-cone distribution A/i for the basic set of 
parameters and three 9 C values. 

Then, for Y > \/2 the 02 mode is the dominant, and for 
Y > 1.8 the X3 mode may dominate (if it is unstable for 
a particular set of parameters). 

We note an interesting feature specific for our electron 
distribution: the amplification of the third harmonics of 
both extraordinary and ordinary waves is ineffective for 
hard electron spectra (e.g., 7 = 6, left column) for most 
of the Y range, while it is effective for softer spectra (e.g., 
7 = 8, right column). The reason for the instability sup- 
pression for the hard spectrum is the destructive contri- 
bution of (true absorption by) high-energy electrons. For 
softer electron distributions the role of high-energy elec- 
trons is less important allowing the instability at the third 
harmonics. Note that for the gaussian angular distribu- 
tion used in (Fleishman & Yastrebov 1994a, 1994b), the 
instability suppression by high-energy tail of power-law 
distribution is important even at the fundamental and sec- 
ond harmonics, since the gaussian distribution provides 
weaker angular gradients than considered here gaussian 
loss-cone distribution (6). 

6.3. Effect of angular distribution. 

The dependences of the maximized (over the emission 
frequency and angle) spatial growth rates on 9 C displays 
bell-shaped curves for any mode. The most effective am- 
plification of Ol waves occurs around 9 C sa 45°, while 
around 9 C sw 60° for 02 waves and even higher values for 
X2 waves. No qualitative difference is seen for various 
electron spectral indices 7. 

Of particular importance is the dependence on the 
width of the distribution Afi, figure 8. Evidently, all 
the growth rates increase as A/i decreases because this 
corresponds to stronger angular gradients, and, thus, to 
stronger instability. However, the growth rates decrease as 
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Fig. 9. Spatial growth rates for X2-, Ol-, and 02- modes vs 
minimal momentum of fast electrons po for the basic set of 
parameters and three 8 C values. 

A/i increases first gradually and then rapidly. This rapid 
decrease occurs when the instability approaches its thresh- 
old provided by negative contribution of the momentum 
term (df/dp = —7/) in the expression for the GS absorp- 
tion coefficient. When A/i exceeds 0.1 — 0.3 (depending on 
other parameters) the instability disappears. We empha- 
size that this happens when the pitch-angle distribution 
of fast electrons remains clearly anisotropic. Thus, such 
anisotropic distributions can generate standard GS emis- 
sion without coherent ECM emission. 

The detailed shape of the curves depends on other pa- 
rameters, like Y and 9 C . For X2 and 02 waves, large 9 C 
values are preferable, while for Ol waves - lower 9 C values. 
However, for all considered cases, 9 C = 60° has the largest 
threshold A/i value. 

6.4- Effect of distribution over momentum. 

The growth rates of all considered wave-modes display 
rather weak dependence on pi> r for the selected set of ba- 
sic parameters, so we do not present the corresponding 
plots. We should note, however, that this property of 
the ECM emission is model-dependent and valid if the in- 
stability operates noticeably above the threshold. If the 
instability is very close to its threshold, the behavior of 
the momentum distribution in the high-energy range may 
be important. In particular, this is the case of gaussian 
pitch- angle distribution (Fleishman & Yastrebov 1994a, 
1994b). 

The dependence of the growth rates on the minimal mo- 
mentum po is shown in figure 9. Shape of the curves is 
clearly different for various harmonics and wave- modes. 
X2 waves display decreasing curves with increasing j»q, 
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Fig. 10. Spatial growth rates for X2-, X3-, Ol-, 02-, and 
03- modes vs spectral index of fast electron distribution 7 for 
the basic set of parameters and three 9 C values. 

thus, low-energy electrons are preferable to produce these 
waves. Wide loss-cones (8 C = 60°, 75°) provide larger 
growth rates than narrower loss-cone (9 C = 45°). 

Ordinary waves are generated most effectively if the 
low-energy cut-off is about 20 keV for 02-mode and about 
10 keV for 02 mode, although the corresponding peaks 
are rather broad. Narrow loss-cone (8 C = 45°) is the most 
efficient to produce Ol waves and wider one (8 C = 60°) 
is preferable to produce 02 waves. The wide loss-cone 
of 9 C = 75° is less effective for the both harmonics of the 
ordinary waves. 

Considering together the plots in figure 9, we can con- 
clude that different fractions of fast electrons are respon- 
sible for generation of these different wave-modes. Thus, 
quasilincar relaxation against one (the dominant) mode 
should not necessarily provide saturation of all other wave- 
modes in agreement with numerical quasilinear simula- 
tions (Fleishman & Arzncr 2000). 

The dependences of the growth rates on the electron 
spectral index 7 are of particular importance, since they 
can be compared with observations of accompanying con- 
tinuum (microwave and HXR) emissions. All the curves 
shown in figure 10 are bell-shaped; peak of each curve 
corresponds to a favorable 7 value specific for each wave- 
mode. The decrease of the growth rates for harder spec- 
tra (smaller 7) is provided by negative contribution of 



high-energy electrons at large harmonics. The decrease 
for softer spectra (larger 7) is related to the increase of 
the instability threshold that is proportional to 7 for the 
power-law momentum distribution of the electrons. 

The favorable 7 value depends on the harmonic num- 
ber and type of the wave. X2 waves have the peak around 
7 = 6 (which corresponds to the synchrotron spectrum 
with the index 1.5) for a broad range of the involved pa- 
rameters. A bit lower values are favorable for 02 waves. 
However, Ol waves display the peak at rather small values 
7 = 4-5. 

It is remarkable that the favorable 7 values for the third 
harmonics are significantly larger than for the fundamen- 
tal and second harmonics: it is about 8 (synchrotron spec- 
tral index 2.5) for 03 waves and even larger for X3 waves. 
Such a different behavior of the curves for the different 
waves can in principle be used to specify the contribution 
of various harmonics to the spike radiation. 

7. Discussion 

We conclude the pitch-angle anisotropy to affect signif- 
icantly the radio emission from solar flares. First of all, 
we should mention the strong effect of the anisotropy on 
the intensity and spectral index of the emission in the op- 
tically thin region, which relates primarily to the energy- 
dependent directivity (1) of the synchrotron radiation by 
a single electron. 

Currently, it is discovered with the NoRH data that 
microwave emission from the foot-point parts of ex- 
tended flaring loops displays systematically softer spectra 
(Aa ~ 1-1.5 in the domain 17 GHz to 34 GHz) than from 
the loop-top part (Yokoyama et al. 2002; Mclnikov 2005). 
Evidently, the pitch-angle anisotropy can make an impor- 
tant contribution to the observed spectral index varia- 
tions. Indeed, if we take into account the convergence of 
the magnetic field lines towards the foot-points, we should 
expect a more anisotropic electron distribution (due to a 
larger loss-cone) in the foot-point source than in the loop- 
top source. Furthermore, for the disc flares the foot-point 
source is observed at a QP direction, while the loop-top 
source is observed at a QT direction, which results in sys- 
tematically softer spectra of GS radiation from the foot 
points compared with the loop-top. Thus, the pitch-angle 
anisotropy probably plays a key role in the observed vari- 
ations of the spectral index of the microwave radiation 
along the loop. Moreover, the pitch-angle anisotropy has 
strong effect on polarization properties of GS emission, 
which is widely observed with NoRH. 

Coherent GS emission (ECM) is also highly relevant for 
understanding the current radio observations. The avail- 
able data about the radio spikes are pretty consistent with 
the idea that the source of spike cluster is a loop filled by 
fast electrons and relatively tenuous background plasma 
(Fleishman et al. 2003). The trapped fast electrons (with 
a power-law energy spectrum and a loss-cone angular dis- 
tribution) produce continuum GS emission through the 
entire loop, while each single spike is generated in a local 
source inside this loop by ECM mechanism when some fa- 
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vorable conditions are fulfilled. Most probably, the local 
source is formed when the local anisotropy is increased 
compared with the averaged one to produce ECM emis- 
sion (Fleishman & Mclnikov 1998). The assumed fluctu- 
ations of the pitch-angle distribution of fast electrons can 
be produced by the magnetic turbulence (Fleishman & 
Melnikov 1998; Barta k Karlicky 2001). Accordingly, if 
the averaged angular distribution is stable against ECM 
generation, the number of spikes in an event should in- 
crease as the intensity of the magnetic turbulence in- 
creases (and other equal conditions). Complementary, the 
closer the averaged angular distribution to the instability 
threshold, the larger the amount of the local sources (and, 
respectively, the number of spikes) might be formed by the 
same magnetic turbulence. 

The scale of a single spike source is estimated (indi- 
rectly) as ~ 200 km (Benz 1985); a wave should expe- 
rience many e-folding amplifications over this length to 
produce observable spike with large brightness tempera- 
ture, so this requirement reads k>5 - 10 _8 cm _1 , where k 
is the spatial amplification coefficient. Wc conclude (from 
the numbers in figure 7) that the required amplification 
can easily be provided for the fundamental and second 
harmonics of both X- and O- modes. Moreover, the am- 
plification of X3 mode is also possible if fast particles are 
numerous enough, N e /N ,c 10 -2 , while the amplification 
of 03 mode is less probable. 

The dependence of the growth rate on the fast electron 
distribution is of primary importance since it can be com- 
pared with the observations. Indeed, the decrease of the 
growth rates with the 7 increase corresponds to smaller 
efficiency of the ECM emission for softer spectra of fast 
electrons. The correlation of this type is actually observed 
when spectra of microwave or HXR emission simultaneous 
to spike bursts are analyzed: the harder the fast electron 
spectrum the stronger the averaged spike flux (Fleishman 
et al. 2003). The range of microwave spectral indices is 
typically a = 1 — 4. Thus, we can conclude that most of 
the spikes are generated at a harmonics not larger than 
second. ECM emission at the third harmonics (see Figure 
10) is the most efficient if 7 r* 8 (a ss 2.5), therefore, the 
softer microwave spectra should be preferable for the third 
harmonics generation in the range a < 2.5, which has not 
been observed for spikes yet. 

Then, the dependences on the angular gradient value 
(A/i) display clearly that the instability disappears (as 
Afj, increases) when the pitch-angle distribution still re- 
mains highly anisotropic (A/x = 0.2 — 0.3). Accordingly, 
the overall distribution of fast electrons accumulated at 
the loop may be rather anisotropic for a long time provid- 
ing noticeably softer GS spectra (for a QP viewing angle) 
than the isotropic distribution, which is observed indeed 
(Fleishman et al. 2003). 
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